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We investigate the stability of the collinear ferromagnetic state in kinetic exchange models for 
(III,Mn) V semiconductors with randomly distributed Mn ions . Our results suggest that noncollinear 
ferromagnetism is commom to these semiconductor systems. The instability of the collinear state 
is due to long-ranged fluctuations involving a large fraction of the localized magnetic moments. 
We address conditions that favor the occurrence of noncollinear groundstates and discuss unusual 
behavior that we predict for the temperature and field dependence of its saturation magnetization. 



In systems with local moments, interactions that fa- 
vor parallel spin alignment (ferromagnetic interactions) 
yield ferromagnetic ground states. Interactions that fa- 
vor antiparallel spin alignment (antiferromagnetic inter- 
actions), on the other hand, yield antiferromagnetism 
only in the special case of short-range interactions on 
bipartite lattices. More generally, longer range antifer- 
romagnetic interactions often lead to complex order or, 
when spatial disorder is also present, to spin-glass ground 
states Q. 

In this context, the occurrence of ferromagnetism 0] in 
III-V semiconductors with Mn substituted on a relatively 
small fraction of the cation sites is unusual, since Mn 
provides S — 5/2 local moments in this system and the 
distribution of Mn ions among the cation sites is thought 
to be random. A ground state with parallel alignment 
can be understood in a system with randomly located 
local moments only if interactions are dominantly ferro- 
magnetic, a property not normally associated with the 
carrier-mediated interactions thought |3|] to be responsi- 
ble for ferromagnetism in these systems. Indeed there is 
ample evidence that randomness can play an important 
role in these ferromagnets, since transport and magnetic 
properties depend strongly on the MBE growth condi- 
tions used for materials synthesis, with reproducible re- 
sults achieved only if growth protocols are followed care- 
fully §]. 

An important step toward the systematic analysis of 
such effects was taken very recently by Potashnik et al. H] 
who studied the effect of annealing procedures on mag- 
netic and transport properties for samples with a Mn 
fraction close to that at which the maximum Curie tem- 
peratures (T c ) occurs, x ~ 0.05. They find that T c , the 
ground state saturation magnetization M(T = 0), and 
the shape of the M(T) curve all depend on the tem- 
perature and duration of annealing. Since the Mn con- 
tent of the system is not believed to be changed by an- 
nealing, changes in the T = magnetization can be ex- 
plained only if Mn moments are not always, or perhaps 
never, completely aligned in the ground state. More- 
over, drastic changes in the shape of the M(T) curve are 



evidence of changes in the magnetic excitation spectrum. 
Both changes most likely reflect the sensitivity of carrier- 
mediated interactions between the Mn local moments to 
the defect distribution that is varied by annealing. 

In this Letter we study the stability of states with 
parallel Mn local moments using the kinetic- exchange 
model of carrier-induced ferromagnetism. We find that 
fully aligned ferromagnetic ground states do not occur 
when the band electrons wavefunctions are strongly con- 
centrated around Mn ion locations, and propose that 
this circumstance must be avoided to achieve high fer- 
romagnetic transition temperatures. On the other hand 
we find, consistent with the observations of Potashnik 
et al., that noncollinear ferromagnetism is common in 
(III,Mn)V semiconductors. On the basis of calculations 
of external-field dependent magnetic excitation spectra, 
we discuss unusual temperature and field dependence of 
the magnetization that we predict will be characteristic 
of non collinear states in these ferromagnets. 

Ferromagnetism in (III,Mn)V semiconductors has been 
addressed using first-principles local-spin-density ap- 
proximation calculations |5j] , phenomenological impurity- 
band models || and, most-commonly [|||7]-§j], using a 
semi-phenomenological approach that sidesteps atomic 
scale physics by using an envelope function description 
of the bands. The last approach is strongly supported 
by the large literature on the closely related (II,Mn)VI 
semiconductors (|10| and we follow it here. The model 
Hamiltonian we study, 
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describes noninteracting carriers in a parabolic band 
characterized by an effective mass m, whose spin den- 
sity s(f) is coupled to ionized Mn acceptor spins Si at 
locations Ri by an spatially extended antiferromagnetic 
exchange coupling J(r), which we take to be of the form 
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where J p d parametrizes the strength and ao the range of 
the interaction. Similar models have been studied some 
time ago |ll| in connection with classical spin glass sys- 
tems; in that case however the density of magnetic im- 
purities N is taken to be much smaller than the density 
of charge carriers n, n/N ^$> 1 pi] . (III,Mn)V semicon- 
ductors are in the opposite, n/N <C 1 because of strong 
compensation of Mn acceptors by antisite defects N. The 
model (uh is unrealistic in its use of a single parabolic 
band to represent the semiconductor valence bands, in 
its neglect of Coulomb interactions between holes in the 
valence band and other holes and Mn acceptors, and in 
the absence of any representation of the defects (other 
than the Mn ions) whose distributions presumably vary 
during the annealing process. We return to these limi- 
tations of the model later. In addition, the range of the 
exchange interaction in this model accounts |1^] for the 
atomic scale physics only crudely. We believe, neverthe- 
less, that this relatively simple model captures much of 
the physics associated with the approach of and the oc- 
currence of reduced saturation magnetization states for 
non-optimal defect configurations. 

Recently a theory of magnetic fluctuations around a 
fully aligned state of magnetic ions was developed us- 
ing the virtual crystal approximation J7|. In contrast to 
that work we do not adopt this simplifying approxima- 
tion in which the Mn magnetic distribution is replaced by 
a continuum, retaining individual spins of length S = 5/2 
placed at arbitrary locations instead. The virtual crystal 
approximation, which has a long and successful history 
in analyzing the properties of (II,Mn)VI ferromagnets, 
completely eliminates disorder and, as we discuss further 
below, noncollinear ferromagnetism. 

We first address the stability of the collinear state in 
which all Mn moments have parallel orientation taken 
to be in the z direction, parameterizing changes in the 
orientation of spin / by a complex number, zi = (Sf + 
iSj)/v2S. The zeroth-order band electron Hamiltonian 
is given by Eq. p) with Si = Sz. For small changes in 
spin orientation the change in the band electron Hamil- 
tonian is 
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The contribution to the carrier groundstate energy of first 
order in the zi vanishes. At second order one obtains 



e(2) = E zi M u z J 
i, j 

where Mu = Lu + Ku and 
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Here (s(r)) is the band-electron spin density in the 
collinear state, n F is the Fermi function at zero tem- 
perature In these equations we have defined 
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with ip a(J {f) being the spin a component of the the car- 
rier wave function with energy e a = £ Q + /i, where [i is 
the Fermi energy. All eigenvalues of M must be posi- 
tive for the collinear state to be stable. The diagonal 
contribution to M from the positive definite matrix (|5|) 
gives the energy change associated with reorienting each 
localized spin in the mean-field created by the polarized 
band electrons, while the contributions from the matrix 
(g) represents the energy reduction due to the response 
of band electrons to reorientation of the localized spins. 
When the Mn spins are treated quantum mechanically 
using a Holstein-Primakoff boson representation [1J] and 
coherent state path integrals the zj are boson coherent 
state labels, eqivalent to replacing zi and zi by bosonic 
creation and annihilation operators aj and a~f . Then 
Eq.( |3|) is the spin- wave Hamiltonian, and the eigenval- 
ues p4| of Mu are the system's spin-wave energies. 
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FIG. 1. Upper panel: Disorder-averaged density of states 
for spin-wave excitations. These results were obtained with 
a simulation cube of volume V — L — 400nm with a Mn 
density of TV = 1.0nm~ 3 , a carrier density n — 0.15nm -3 , 
and m = 0.5m e ;. The strength of the exchange interaction 
between ions and carriers is J p d = 0.05eVnm with a spa- 
tial range of ao = 0.40nm. Lower panel: Disorder- averaged 
participation ratios for the same system. In both graphs the 
value at zero energy is enhanced due to the contribution of 
the uniform rotation mode which occurs for any disorder re- 
alization. 

We have evaluated the matrix Mjj and its spectrum by 
solving for the band electron mean-field eigenstates, ap- 
plying periodic boundary conditions to cubic simulation 



cells. The single-particle wavefunctions i[>aa{r) are com- 
puted in a plane-wave basis taking into account wavevec- 
tors q with length up to an appropriate cutoff q c . The 
upper panel of Fig. [l| shows typical results for the eigen- 
value distribution of Lu averaged over a large number of 
randomly chosen Mn distributions. We have established 
by further calculations that the effects of the wavevec- 
tor cutoff on the low-lying excitations have already sat- 
urated for the value of q c used here. We note that for 
any arrangement of the Mn positions Rj the matrix Mjj 
contains a zero eigenvalue corresponding to a uniform ro- 
tation of all spins where all components zi are the same. 
This implies summing over all columns of M gives the 
zero vector, a sum rule that is perfectly fulfilled in our 
numerics independently of q c since this plane wave cutoff 
does not influence the rotational invariance in spin space. 

In the above calculations the Mn positions were cho- 
sen completely at random with a uniform distribution, 
while in a real (III,Mn)V semiconductor the Mn ions are 
located on the fee lattice of cation sites. Our disorder- 
averaged data for the two cases is are indistinquishable, 
however, establishing that atomic-scale correlations in 
Mn locations do not alter our results. Our most impor- 
tant finding, illustrated in Fig. [I] and discussed at length 
below, is that that for a wide range of parameter values 
there are always a few negative eigenvalues, i.e. that the 
collinear ferromagnetic state is not stable. 

To analyze the nature of this instability we have evalu- 
ated participation ratios for these elementary excitations 
which we define by 
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where a\ is the 7-th component of the j'-th normalized 
eigenvector of Mij. pj is an estimate for the fraction 
of Mn sites that have important involvement in the j'-th 
spin wave. For instance, if a vector contains exactly a 
fraction of p nonzero components of equal modulus and 
all others being zero, its participation ratio is p. The 
largest participation ratio of unity is achieved for the 
zero-energy uniform rotation mode where all components 
of the corresponding eigenvector are equal. 

The lower panel of Fig. [j] shows the disorder-averaged 
participation ratio as a function of spin-wave energy for 
the same situation as in the top panel. The property 
that negative energy excitations have large participation 
ratios shows that the instabilities of the collinear state in- 
volve correlated reorientations of many spins, rather than 
lone loosely coupled moments. It is the high energy ex- 
citations, for which the mean-field term dominates, that 
are single-ion in character with participation ratios com- 
parable to 1/400 where 400 is the number of Mn ions 
in these simulations cells. The shape of the spin-wave 
density-of-states is, in the model we have studied, sen- 
sitive to the Mn density N, the carrier density n, and 



to the Hamiltonian parameters m, J p d, ao- Situations in 
which the collinear ferromagnetic state is, for certain dis- 
order realizations, stable can be approached most simply 
(for technical reasons) by letting ao be larger, but also 
we believe for larger n/N. For the value of n/N illus- 
trated in Fig. [I], which corresponds to a carrier density 
somewhat smaller than that measured in the highest T c 
samples, negative eigenvalues occur for nearly any Mn 
ion distribution. We note that negative eigenvalues in- 
crease in number as the wavevector cutoff is increased 
toward its converged value. These results suggest that 
noncollinear ferromagnetic states are common, and that 
they are sensitive to the distribution of Mn ions and de- 
fects, especially those that trap carriers. We find that the 
collinear ferromagnetic state tends to become unstable as 
mean-field band eigenfunctions become more strongly lo- 
calized around Mn ion sites. 

The effect of a weak external magnetic field on the 
eigenvalue spectrum of Mu is particularly simple. The 
field couples to the local moment through its Lande g- 
factor, adding 2fijsH to the energy of a spin-wave exci- 
tation for S = 5/2, and to the spin and orbital degrees of 
freedom of the band electrons. The orbital coupling leads 
to Landau levels that do not play an important role at 
weak fields in these highly disordered samples. Zeeman 
coupling is also unimportant, because it gives a contribu- 
tion to the spin-splitting that is negligible compared to 
the mean-field exchange splitting A = J p dNS ~ O.leV. 

Thus, provided that the modulus of the smallest eigen- 
value (i.e. the onset of the eigenvalue distribution of 
Mu at negative energies) is small compared to A (as 
it is the case in data shown above), this modulus is the 
value of gLfJ-BH necessary to force full spin alignment of 
a noncollinear state. The experiments of Potashnik et 
at demonstrate that the maximum value of M(T = 0) 
is achieved over a certain range of annealing histories. 
We associate this maximum value with the fully aligned 
collinear Mn configuration state; indeed the maximum 
moment per Mn is consistent with full alignment par- 
tially compensated by band electrons. In this circum- 
stance, the virtual crystal approximation used in earlier 
spin-wave calculations which neglected disorder (jj can 
be accurate. For other annealing histories, M(T — 0) is 
reduced corresponding, we propose, to noncollinear or- 
der of the Mn spins. Calculations like those described 
above show that these spins gradually align as an exter- 
nal field is added to the Hamiltonian. Since full align- 
ment is achieved at a finite external field, we expect full 
alignment to be indicated experimentally by a kink in the 
M(T = 0,H) curve. At this point, we predict that the 
system will still have gapless excitations and power-law 
temperature dependence of the magnetization, in sharp 
contrast to the gapful excitations and exponentially sup- 
pressed temperature dependence of conventional ferro- 
magnets in an external magnetic field. At the same 
time, spin-resonance experiments will see a gapped spin- 



wave spectrum since they couple only to the uniform- 
rotation mode which has finite energy g^nsB over the 
noncollinear groundstate. 

To address the noncollinear state, we generalize our 
formalism by expanding around and defining Holstein- 
Primakoff bosons pl| with respect to a general non- 
collinear orientation configuration, of each Mn spin, Si — 
Sttj = S (smtf i cos ipi, s'mfli sin (pi, cos di). If this con- 
figuration is not an extremum, there will be an energy 
correction that is linear in the Holstcin-Primakoff vari- 
ables: 
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where e Vl — (cos<£>/,sin<^j,0) and e z — (0,0,1). 
The contributions to the carrier ground state energy 
quadratic in the Holstein-Primakoff variables can also be 
obtained via perturbative calculations similar to those 
for the the collinear case. The components gj represent 
the gradient of the energy with respect to deviations, pa- 
rameterized by the zj, from the reference orientations. 
A given orientation of Mn spins is stationary with re- 
spect to fluctuations if all complex coefficients gi vanish. 
This is the case if and only if J d 3 rJ(r — Rj)(s(f)) is 
parallel with the direction Clj of the local ion spin. (The 
collinear ferromagnetic state is always stationary but not 
necessarily stable.) 

We have employed the energy gradient expression (||) 
in a numerical steepest descent procedure to search for 
true energy minima. Briefly, our results are as follows. In 
cases where the energy minimum found by this method 
is close to that of the collinear state (with M(T = 0) 
about 90 or more percent of the maximum value), this 
minimum appears to be unique for each disorder realiza- 
tion. We can therefore be confident that we have located 
the absolute ground state of our model. In situations 
where the magnetisation is reduced more substantially, 
however, by about 20 percent or more typically, we con- 
verge to different energy minima from different starting 
points. In these case the model has substantial spin-glass 
character with a complex energy landscape. For the sys- 
tem shown in Fig. |l|, for instance, magnetisation values 
at local energy minima are typically 30 to 40 percent of 
the collinear state value. 
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